Microscopic model for the logarithmic size effect on 
the Curie point in Barabasi- Albert networks 

Agata Aleksiejuk-Fronczak 
February 6, 2008 



Institute for Theoretical Physics, Cologne University, 
D-50923 Koln, Euroland 

and 

Faculty of Physics, Warsaw University of Technology, 
Koszykowa 75, PL-00-662 Warsaw, Poland [} 

Abstract 

We found that numbers of fully connected clusters in Barabasi- 
Albert (BA) networks follow the exponential distribution with the 
characteristic exponent k — 2/m. The critical temperature for the 
Ising model on the BA network is determined by the critical temper- 
ature of the largest fully connected cluster within the network. The 
result explains the logarithmic dependence of the critical temperature 
on the size of the network N. 
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During the last few years studies of random, evolving networks have be- 
come a very popular research domain among physicists. A lot of efforts were 
put into investigation of such systems in order to recognize their structure 
and to analyze emerging complex dynamics. The Barabasi-Albert (BA) 
network [|J being the subject of this paper is probably the most studied 
model that realizes many properties of real weblike systems. The most im- 
portant feature of the model is the power-law degree distribution that is 
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also a characteristic of many real networks. It was shown that this feature 
is closely related to processes governing system evolution. Two important 
ingredients of the BA model are: continuous network growth and preferen- 
tial attachment. The network starts to grow from an initial cluster of m 
fully connected sites. Each new node that is added to the network creates 
m links that connect it to previously added nodes. The preferential attach- 
ment means that the probability of a new link to end up in a vertex i is 
proportional to connectivity of this vertex k{. 

The goal of the investigation run here was to understand the logarithmic 
dependence of the Curie temperature in ferromagnetic phase transition in 
BA network § on the system size. A mean field theory for this behaviour 
has been recently found by Bianconi || and independently developed by 
Dorogovtsev at al. Q|. Here we provide another explanation for the fact and 
associate it with the microscopic, internal structure of the BA networks. 

We investigate the probability to occur a fully connected cluster of size s 
in BA networks. A fully connected s-cluster is defined as a group of s points, 
such that every point of the group has a connection to each of the (s — 1) 
remaining points within this group. In order to get reliable numbers n s of 
fully connected s-clusters all s-subsets of the whole set of network vertices 
should be examined. It gives ) possibilities for each cluster of size s and 
altogether 2^ combinations for all cluster sizes. Such calculations would 
involve unreasonable amount of computer time. To simplify the task we 
noticed that each point i that has k{ nearest neighbors may be involved in 
(s-i) f un y connected clusters of size s. Because of computer time limitations 
we were forced to restrict our studies to nodes with degree h L < k max where 
kmax = 20 or 30 and small lattices ./V = 28, 38, 53 (Fig.[l]). All simulations 
were done for BA networks with m = 3. Results presented at Figjl] are 
exact in this sense that all possible combinations of fully connected clusters 
up to size s = k max + 1 having at least one node i of degree h L < k max were 
checked. Numerical simulations carried out for networks of size ./V show that 
the cluster numbers n s follow the exponential distribution: 

n s = ANe~ K * s . (1) 

At the Fig.|| numbers of such s-clusters per network node n s /N are fitted 
to @j giving 

hx{A) = 3.20 ± 0.20 (2) 
k = 0.70 ± 0.05. 

The critical temperature for the Ising model on the fully connected lattice 
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that consists of s sites follows the well-known mean-field result 

T c = {s- 1), (3) 

where (s — 1) is the mean number of nearest neighbors (temperature is 
measured in units of coupling constant over Boltzman constant J/k). Let 
us assume that the critical temperature for the largest fully connected cluster 
within the BA network is given by (]3|). Since fully connected clusters in BA 
model are not isolates lattices but they are involved in the whole network 
structure the equation (^) is valid only for the largest cluster. Putting n s = 1 
and using (|l|) we get the temperature as: 

ln(A) ln(Af) 
T c = — - 1 + — (4) 

K K 

Taking into account the result (||) for the BA model with m = 3 we get: 

T c = Bln(N) + C, (5) 

where B = 1.4 ± 0.1 and C = 3.6 ±0.1. Looking at the result obtained by 
Dorogovtsev at al. Q for the whole BA network: 

T c = m\n{N)/2 (6) 

and putting m = 3 into (^) the critical temperature of the whole network 
gains the form: 

T c = 1.51n(iV), (7) 

that for large N fairly good agrees with (|5|) and (Q) giving k = 2/m. It fol- 
lows that the Curie temperature of the whole BA network is closely related to 
the Curie temperature of the largest fully connected cluster. To destroy the 
magnetisation of the network one needs to destroy ferromagnetic coupling 
within this cluster. Our studies are consistent with the general conclusion of 
the paper Q that the most connected vertices in random graphs with arbi- 
trary degree distributions induce strong ferromagnetic correlations in their 
close neighborhoods at high temperatures. Such vertices probably partici- 
pate in the fully connected cluster making the ordering phenomena of the 
whole network dependent on ordering of the largest fully connected cluster. 

In summary, we presented a possible microscopic explanation for loga- 
rithmic size effect of the Curie temperature in Barabasi-Albert networks. 
We found that the critical temperature of the whole network is determined 
by the critical temperature of the largest fully connected cluster within the 
network. 
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Figure 1: Distributions of fully connected clusters, k max = 30 (a) or 20 (b) 
(averaged over 50 networks) 
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